Entangling power of baker's map: Role of symmetries 
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The quantum baker map possesses two symmetries: a canonical "spatial" symmetry, and a time- 
reversal symmetry. We show that, even when these features are taken into account, the asymptotic 
entangling power of the baker's map does not always agree with the predictions of random matrix 
theory. We have verified that the dimension of the Hilbert space is the crucial parameter which 
determines whether the entangling properties of the baker are universal or not. For power-of-two 
dimensions, i.e., qubit systems, an anomalous entangling power is observed; otherwise the behavior 
of the baker is consistent with random matrix theories. We also derive a general formula that relates 
the asymptotic entangling power of an arbitrary unitary with properties of its reduced eigenvectors. 

PACS numbers: 03.65.Ud, 03.65.Yz, 05.45.Mt 



INTRODUCTION 



The baker's map was quantized in 1987 by Balazs and 
Voros PI and soon became a very useful toy model for 
investigating quantum-classical correspondence issues in 
closed chaotic systems, like the scarring phenomenon, 
Gutwiller trace formula and the long-time validity of 
semiclassical approximations (see, e.g., 0, 0, H, H, IE 13 ) ■ 

Later on, the quantum baker appeared in a variety of 
problems of Quantum Information, Quantum Computa- 
tion and Quantum Open Systems. Schack noted that 
the quantum baker could be efficiently realized in terms 
of quantum gates Q. A three qubit Nuclear Magnetic 
Resonance experiment was proposed p and then imple- 
mented (with some simplifications) 1(J|. On the theoreti- 
cal side, Schack and Caves 11] showed that the quantum 
baker of Balazs and Voros' can be seen as a shift on a 
string of quantum bits -in full analogy with the classical 
case- and exhibited a family of alternative quantizations. 
This fainily of bakers was the subject of several studies 
[Tl. [13I IT^ j . Decoherent variants of the baker map have 
been constructed by including mechanisms of dissipation 
and/or diffusion [lElIEIIl- 

The ability of the baker family to generate entangle- 
ment was studied by Scott and Caves 14]. They con- 
cluded that "the quantum baker's maps are, in general, 
good at creating multipartite entanglement amongst the 
qubits. It was found however, that some quantum baker's 
maps can, on average, entangle better than others, and 
that all quantum baker's maps fall somewhat short of 
generating the levels of entanglement expected in ran- 
dom states. This might be related to the fact that spatial 
symmetries in the baker's map allow deviations from the 
predictions of random matrix theory." 



The purpose of the present paper is to demonstrate 
that the spatial symmetry is not to blame for the reduced 
entangling power of the baker. Two numerical comple- 
mentary proofs will be presented. First we check that if 
the symmetry is removed from the baker (by block diago- 
nalization), the resulting desymmetrized bakers produce 
the same levels of entanglement as the original one. Sec- 
ond we show that the entangling power of an ensemble 
of "spatially-symmetric" unitary operators is not signif- 
icantly different from that of the CUE ensemble [ig of 
random operators, i.e., imposing symmetry does not re- 
duce the entangling power. We complete the analysis by 
verifying the dimension of the Hilbert space is indeed the 
crucial parameter which determines whether the entan- 
gling properties of the baker are universal or not. For 
qubit systems, i.e., power-of-two dimensions, an anoma- 
lous entangling power is observed; otherwise the behavior 
of the baker is consistent with random matrix theories. 

The background of this contribution is the wider prob- 
lem of understanding what makes a unitary operation a 
good entangler. There is not a definitive answer to this 
question yet, but some partial results have been obtained 
recently (see, e.g., pj, 2_0, 21, 22J). 

Section ^presents the measure we use for quantifying 
the entangling power of a unitary operation. The quan- 
tum maps to be considered are introduced in section ITTll 
and their symmetries analyzed. Sections IIVI and con- 
tain the numerical analysis of the influence of the sym- 
metries on the entangling power of the baker. The gen- 
eral relation between asymptotic entangling power and 
eigenvectors is discussed in IVII Concluding remarks are 
presented in section IVlTI 



II. ENTANGLING POWER 
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We will only consider the case of bipartite entangle- 
ment of pure states. 

A full system, with Hilbert space of dimension d = 
cIa ^ ds, is partitioned into two subsystems, A and B, 
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with dimensions (Ia and d^, respectively, and such that 
the fuU space Ti has the structure of a tensor product, 
H = Ti.A ^ Ti-B- A usual definition of entanghng power, 
ep{U), of a unitary operator U defined on Ti. relies on a 
vector entanglement measure and on a suitable aver- 
age over initial states p^ : 



ep(C/) = (£;(C/|V'A>®|^B)))|^,),l^, 



(1) 



Thus, ep{U) says how much entanglement U produces, 
on average, when acting on a set of non-entangled states. 
Here we take the measure E to be the linear entropy of 
the reduced density matrix: Let \ip) be a (pure) separa- 
ble state of the full system, corresponding to the density 
matrix p = \ip){ip\. In general, after application of U, the 
new density matrix p' = UpW will not correspond to a 
separable state any more. This will manifest in a positive 
linear entropy of the reduced density matrices 



Sr 



1-tr {pUf ^l-trip'sf >0 



(2) 



where p'^ = tr bp' and p'g — tr ap' |2J|. 

For the average over product states, indicated by 
{■ ■ ■)\i'A),\4>B) CDj we choose the unitarily invariant 
measures in both TCa and Hb [2^ l25l |. That is, the 
components of IV'a) and \iPb) have the same distribution 
as the columns of CUE matrices of dimension x cIa 
and (Ib x d^, respectively p^ . 

Among various possible definitions of the entangling 
strength of a unitary H El IM El , we chose Eq. ((TJ 
mainly for the purpose of comparing our results with 
those in Ref. where that definition was adopted. 

The use of Sl instead of the more natural von Neumann 
entropy not only does not lead to qualitative differences 
[T^ but has the essential advantage of allowing analytical 
calculations, which will be important for understanding 
and extending our results. 



III. QUANTUM MAPS, SYMMETRIES 

Following Schack and Caves we write the unitary 
operator for the quantum baker on N qubits as {d = 2^) 



(3) 



where Gd is the antiperiodic quantum Fourier transform 
on N qubits, I2 is the unit operator for the first qubit, 
and G'^J^^ the inverse (antiperiodic) Fourier transform on 
the remaining N— 1 qubits. We use the standard ordering 
for the computational basis: If \j) is a tensor product of 
individual qubit basis states |ei), with e.; = 0, 1, i.e.. 



\j) = ki) |e2) . . . lew) 
then j is given by the binary expansion 

N 



N-i _ ^ 

= ei . . . Eat 



(4) 



(5) 



< j < d- 1. 

The definition © becomes equivalent to the baker of 
Balazs-Voros and Saraceno when one identifies the com- 
putational basis states \ j) with "position" eigenstates \qj) 
having eigenvalues qj = O.ci ... eat 1. In any case, the ma- 
trix representation of the baker is 



\Bd 



\Gd 



^d/2l 



\G 





d/2l 



(6) 



where ||Gd|| is the inverse Fourier matrix 

The baker map has a time-reversal symmetry, corre- 
sponding to the Fourier transform followed by complex 
conjugation in the computational basis: 



\Gd ^BdGdW* - \\B^ ^1 



(7) 



The use of the antiperiodic Fourier transform makes the 
baker also reflection symmetric, in agreement with its 
classical counterpart 01 ■ That is, if we deflne the reflec- 
tion operator 



then 



Rd\j) = \d-l-j) 



BdRd — RdBd ■ 



(8) 



(9) 



For a qubit system the reflection operator can be factored 
into a tensor product of N single qubit reflections 



Rd\j) ^ R2\ei) ® R2\e2) 



R2\^n) 



(10) 



where R2 is just the negation operator (Pauli-X gate). 
The reflection symmetry of the baker, Eq. ©, can be 
easily proved using the reflection symmetry of the an- 
tiperiodic Fourier transform, RdGd = GdRd, and the 
factorization property Rd = R2® Rd/2- 

Given that the quantum baker is a unitary operator 
with a chaotic classical limit, one may expect that the it- 
erative application of the baker to random non-entangled 
states could produce states with levels of entanglement 
typical of random states. However, Scott and Caves ver- 
ifled that, in spite of being a good entangler, the baker 
generates states that are somewhat less entangled than 
random states. They suggested that this deviation might 
be due to the spatial symmetry Rd, which, as any unitary 
symmetry, is known to produce deviations from standard 
random matrix behavior. 

In the following we implement a simple test to decide if 
the spatial symmetry does really play a significant role in 
the reduction of the entangling power of the baker. Due 
to the symmetry Rd, with eigenvalues ±1, the baker can 
be cast into block-diagonal form: 



At BdA 



where b'^^ are the symmetry-reduced baker maps, and 

At is a unitary mapping of the computational basis to a 
i?(i-symmetrical basis: 

1 



>B 



i-;l + \i){i\®Bi,j> 



(+) 



(11) 



A== {Id + iY Rd/2) , 



(12) 
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with Y the second PauH matrix. 

The maps ^^^2 ^^'^^ '^sU known classical limits: They 
correspond to conservati ve, p iecewise linear versions of 
the Smale horseshoe [s^, l3l( . Instead of working with 
the exact B^^^ , we prefer the simpler approximate ex- 
pressions: 

B^d~^ = Gd(|0)(0|®G-/2 + |l)(l|®Grf/2) 

= Dd, (13) 

RdB^+^Rd - Gd (|0)(0| ® G-]^ ~ |1)(1| ® 

^ D'd- (14) 

The operators on the RHS in the equations above, Dd 
and -D^, are both unitary, time-reversal symmetric, and 
do not have spatial symmetries. They share a same clas- 
sical limit with the reduced bakers, but correspond to a 
slightly different quantization of the classical baker map 
|30|. Instead of Hll|) . D and D' satisfy 

At Brf A = |0)(0|i?rf/2 + \\){\\RdD'^i^Rd , (15) 
where Bd coincides with Schack and Caves' baker Bi [lTI| . 



T « 1 « 1 « 1 « r 



0,9- 




20 40 60 80 100 120 

n 



FIG. 1: (color online) Five initial pure product states were 
chosen randomly according to the CUE measure, and then 
evolved by applying the D-map n times. At each time the 
linear entropy. Shin), is calculated. The system consists of 
eight qubits divided into two subgroups of most/least signifi- 
cant four. 



IV. COMPARING Bi WITH Di 

In order to assess the effect of symmetries in the pro- 
duction of entanglement, it suffices to consider the non- 
symmetric operators Dd or Z?^, and compare with Bd- 

The iterative application of Bd to an initial product 
state typically makes the entropy grow quickly from zero 
to some "equilibrium" value, which depends on the ini- 
tial state. After that, the entropy keeps fluctuating with 
small amplitude around the equilibrium. This is the be- 
havior observed by Scott and Caves for a variety of ways 
of partitioning the qubits [l^ . Our analises of the map 
Dd [Eq. (|13|l ] verify the same qualitative features. As an 
example, we show in Fig.^some plots of entropy vs. time 
for a system of eight qubits split into two groups: the four 
most significant qubits on one side, the remaining least 
significant on the other. 

From now on we focus on the asymptotic regime of 
large times, when the system has already relaxed to equi- 
librium. In this regime we expect the statistical proper- 
ties of entangled states to be described by some random 
matrix model. The simplest ansatz associates evolved 
pure states with random states of the same Hilbert space, 
chosen according to the CUE measure. The average lin- 
ear entropy of these states is given by 



(5 



i/CUE 



{dA-\){dB-l) 



dAdi 



1 



(16) 



analytical expressions for the second and third cumulants 
are also known (see j3| references therein). 

For a quantitative comparison between D-map and 
baker, we considered the same system as before, but this 



time we generated a set of 2 x 10^ data by gathering 
values of SL{n) for 513 <n< 2512 and 10'^ random ini- 
tial states. These data, properly binned, are displayed in 
Fig. [Altogether with analogous data for the baker. It can 
be immediately seen that both maps, baker and D^, pro- 
duce very similar distributions of entropies, both shifted 
to values lower than those of random states. The con- 
clusion of this comparative simulation is that the reflec- 
tion symmetry is not the cause for the states generated 
by the baker being less entangled, in average, that ran- 
dom states (because D-map is non-symmetric and also 
shows a reduced entangling power). However, absence 
of symmetry may be the explanation for a very small, 
though perceptible, increase of entangling power of the 
ZJ-map as compared with the baker (see Fig.[2Jl. (Similar 
influence of symmetry in coupled tops was reported by 
Bandyopadhyay and Lakshminarayan j32].) 

Another simple complementary test, which also shows 
that the influence of the symmetries is very limited, con- 
sists of calculating the entangling power of a random 
map having the same symmetries as the baker, i.e., time- 
reversal and reflection symmetries. 

The problem of introducing symmetries in a random 
matrix model is well known, e.g., in the scattering ap- 
proach to electronic transport through mesoscopic cavi- 
ties [33, 0, : The unitary (scattering matrix) is 
brought to block diagonal form by choosing a basis with 
well defined symmetry, and each block is modeled by a 
circular ensemble. Using this recipe for the baker map 
we arrive at an ensemble of random matrices B with the 
structure 



6 = A 



(|0)(0| ® iy« + |1)(1| ® At ^ (17) 
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FIG. 2: Histograms of asymptotic linear entropies generated 
by the baker map (open circles) and by the map Dd (full 
circles). The line corresponds to the distribution of entropies 
for a large set of CUE random states calculated numerically. 
In all cases the system consists of eight qubits divided into 
two subgroups of most/least significant four. 

where W^*^^^ and VF'^'' are drawn independently from the 
COE ensemble, ap prop riate for unitary maps with time 
reversal symmetry [l^! and A is defined in (|12|) . 

In Fig. Owe calculate the entangling power of the sym- 
metric random maps defined above and compare with 
random maps having no symmetries at all, i.e., the CUE 
ensemble. The figure shows that differences between both 
ensembles are not significant. 

(CUE matrices were generated using the Hurwitz pa- 
rameterization [37L l38l |. COE matrices were obtained 
simply by forming the products VV'^ , with V belong- 
ing to CUE JIJ.) 



V. FROM QUBITS TO ARBITRARY 
DIMENSIONS 

In the search for an alternative explanation for the 
anomalies observed, we recall that baker maps in spaces 
of power-of-two dimensionality are known to exhibit pe- 
culiar properties. For instance, Balazs-Voros and Sano 
observed that bakers of dimension 256 and 1024, respec- 
tively, display spectral statistics quite far from universal 
behavior 2, 39J. Though asymptotic entangling power is 
not a property determined by the eigenvalues, but by the 
eigenvectors (see below), the latter may also be anoma- 
lous for the qubit case. Thus, we now proceed to check if 
the precise value of Hilbert space dimension has a definite 
influence in the entangling properties of the baker. 

Note that Eq. (|SJ| defines a quantum baker for any even 
dimension d. We shall consider, for instance, two sys- 
tems with dimensions d = 238 and d — 162, partitioned 
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FIG. 3: Histograms of linear entropies generated by one ap- 
plication of random maps belonging either to the symmetric 
ensemble of Eq. (1171 (circles), or to the CUE ensemble (line). 
In the first case we applied 1000 symmetric maps to 1000 ran- 
dom non-entangled states. The CUE data set is the same as 
in Fig. 121 and so is the system of qubits. 

as 238 = 14 X 17 and 162 = 9 x 18. The analyses of 
these cases is presented in Fig.0] The histograms clearly 
show that, by avoiding power-of-two dimensions, we re- 
cover universal behavior. These results give additional 
support to the belief that qubit bakers (and D-maps) 
possess hidden symmetries, as happens with cat maps 
and with certain triangles of the hyperbolic plane [i^l . 

The case d = 162 = 2x3'' was chosen to demon- 
strate that a power-of-three factor in Hilbert space di- 
mension is not a source of anomaly for the baker. Pre- 
sumably, power-of-three dimensions, i.e., qutrit systems, 
are anomalous for ternary bakers, like those considered 
by Nonnenmachcr and Zworsky plll4^. 

VI. RELATION WITH EIGENVECTORS 

We have characterized the entangling abilities of a uni- 
tary operator by the distribution of entropies it produces 
when acting iteratively on a set of non-entangled states. 
But such an entropy distribution depends (weakly) on 
the number of iterations, even in the long-time regime. 
For this reason, in the histograms of Figs. |21 and 01 we 
also included data corresponding to different times. 

Then it is natural to consider an average of the entropy 
over both initial states and time. In this way we arrive at 
a quantity that characterizes unambiguously the asymp- 
totic entangling power of an unitary, e^(U), defined by 

1 ^ 

e-(f/)EE hm ^e,([/'=), (18) 

fc=l 

i.e., the time average of the entangling power of Eq. (Q. 
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FIG. 4: Histograms of asymptotic linear entropies generated 
by two non-qubit baker maps with: (full circles) d = 238, 
cLa ~ 14, and ds = 17; (open circles) d = 162, dA = 9, and 
ds = 18. Data were collected by the same procedure used for 
Fig- HI Full lines are the corresponding CUE predictions. 



In the case of the bakers and D-map considered in 
the previous sections can be estimated as the mean 
value of the data set used to generate each histogram 
in Figs. 121 and 0] However, extracting e^{U) from a fi- 
nite data set introduces undesirable statistal errors. This 
could be avoided by implementing the averages in lfTH|) 
analitycally, followed by numerical evaluation of the re- 
sulting expression. We shall see in the following that this 
procedure leads to a relation between the asymptotic en- 
tangling power of an arbitrary unitary U and its reduced 
eigenvectors. Such a relation is not only very interesting 
by itself, but it will give us the possibility of checking the 
numerical simulations of previous sections. 

We start by noting that if the average over initial states 
is removed from Eq. I|18|l one obtains the asymptotic en- 
tropy considered by Demkowicz-Dobrzanski and Kus |22j| : 

^r(i^» = i-"Ei(e.i^>i'tr^ {pif- 



(19) 



which depends of the initial state l^p). In the expression 
above \ei) stands for an eigenvector of U, and is the 
reduced operator 



Pa = trB|e,)(ej| 



(20) 



Of course, Pa is the reduced density matrix obtained from 
the i-th eigenvector, to be denoted simply by p\. 

Formula (|19|l is not of completely general validity. 
However, its derivation makes only the weak assumption 



that eigenvalues exp(iti 
rability relation 



;) do not satisfy the commensu- 







(21) 



except for the trivial cases k = I and m = n, or k = n 
and / = m (this condition is more general that merely 
requiring absence of degeneracies). We have verified nu- 
merically that the eigenvalues of the maps studied in this 
paper, even if anomalous in other sense, are not commen- 
surable. So, we can safely use (|19|l . 

Averaging H19|) over initial random product states we 
arrive at a formula for e^([/) in terms of the eigenvectors 
of U. The derivation is somewhat lengthy but simple. It 
requires averaging products of the type 



(22) 



where Cq are the coefficients that arise from expanding 
the state \ip) in the eigenbasis of U. The coefficients Cq, 
are distributed like the elements in a column of a CUE 
matrix. The average above is one among others calcu- 
lated by Mello some time ago |4^. We omit the details 
and just show the final result: 



1 



7 J2 [^'^'^ (paPa) + (psPs 



(23) 



Here we used the abbreviations d — dAds and d' = (dA- 
l)(c?_B + 1), together with the property 



trA 



{pip'2)^trB{p],p's) , (24) 

showing that (|23|l and (|19|l are indeed invariant with re- 
spect to the swap of subsystems A and B. 

Equation (|23|l is a useful formula which has absorbed 
the averages analitycally; it expresses the asymptotic en- 
tangling power of a unitary as a function of a special 
combination of pairs of eigenvectors. 



tr^ 



(p^aPa) + (pBpi) 



(25) 



i.e., the symmetrized Hilbert-Schmidt scalar product of 
the reduced density matrices. For the maps considered 
here, we have checked that the calculation of using 
either H23|) or by straightforward time and ensemble av- 
erages leads to consistent results; thus we verified the 
correctness of both procedures. 

We remark that is not directly related to the eigen- 
vector entropies 



1 - trA {p\) 



(26) 



even though, bounds relating asymptotic entangling 
power and eigenvector average entropies can be obtained 
by the use of Cauchy-Schwartz inequality (One must 
remember that in some cases eigenvector entanglement 
may give a wrong estimatimation of [2^l43 |.') 
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VII. CONCLUDING REMARKS 

We demonstrated that the deviations from universal 
behavior reported by Scott and Caves are not due to 
spatial or time-reversal symmetries. Instead, the anoma- 
lous entangling power of the qubit baker originates from 
specificities associated to the dimension of the Hilbert 
space being a power of two. When other dimensions 
are considered, a behavior consistent with random ma- 
trix theory is recovered. Presumably qubit bakers pos- 
sess symmetries of number theoretic origin, i.e., with no 
classical analogues ("pseudo-symmetries" |4^). All the 
members of the Schack-Caves family (of which the 



baker considered here is a special case) suffer, to differ- 
ent extent, from a reduction of the entangling power |14| . 
It is tempting to speculate that those differences may 
be related to each member having a different number of 
pseudo-symmetries . 
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